The narrow footprint of ancient balancing selection revealed by heterokaryon incompatibility genes in Aspergillus fumigatus.

In fungi, fusion between individuals leads to localized cell death, a phenomenon termed heterokaryon incompatibility. Generally, the genes responsible for this incompatibility are observed to be under balancing selection resulting from negative frequency-dependent selection. Here, we assess this phenomenon in Aspergillus fumigatus, a human pathogenic fungus with a very low level of linkage disequilibrium as well as an extremely high crossover rate. Using complementation of auxotrophic mutations as an assay for hyphal compatibility, we screened sexual progeny for compatibility to identify genes involved in this process, called het genes. In total, 5/148 (3.4%) offspring were compatible with a parent and 166/2142 (7.7%) sibling pairs were compatible, consistent with several segregating incompatibility loci. Genetic mapping identified five loci, four of which could be fine mapped to individual genes, of which we tested three through heterologous expression, confirming their causal relationship. Consistent with long-term balancing selection, trans-species polymorphisms were apparent across several sister species, as well as equal allele frequencies within A. fumigatus. Surprisingly, a sliding window genome-wide population-level analysis of an independent dataset did not show increased Tajima's D near these loci, in contrast to what is often found surrounding loci under balancing selection. Using available de novo assemblies, we show that these balanced polymorphisms are restricted to several hundred base pairs flanking the coding sequence. In addition to identifying the first het genes in an Aspergillus species, this work highlights the interaction of long-term balancing selection with rapid linkage disequilibrium decay.

A life-history allele of large effect shortens developmental time in a wild insect population.

Developmental time is a key life-history trait with large effects on Darwinian fitness. In many insects, developmental time is currently under strong selection to minimize ecological mismatches in seasonal timing induced by climate change. The genetic basis of responses to such selection, however, is poorly understood. To address this problem, we set up a long-term evolve-and-resequence experiment in the beetle Tribolium castaneum and selected replicate, outbred populations for fast or slow embryonic development. The response to this selection was substantial and embryonic developmental timing of the selection lines started to diverge during dorsal closure. Pooled whole-genome resequencing, gene expression analysis and an RNAi screen pinpoint a 222 bp deletion containing binding sites for Broad and Tramtrack upstream of the ecdysone degrading enzyme Cyp18a1 as a main target of selection. Using CRISPR/Cas9 to reconstruct this allele in the homogenous genetic background of a laboratory strain, we unravel how this single deletion advances the embryonic ecdysone peak inducing dorsal closure and show that this allele accelerates larval development but causes a trade-off with fecundity. Our study uncovers a life-history allele of large effect and reveals the evolvability of developmental time in a natural insect population.

The human fungal pathogen Aspergillus fumigatus can produce the highest known number of meiotic crossovers.

Sexual reproduction involving meiosis is essential in most eukaryotes. This produces offspring with novel genotypes, both by segregation of parental chromosomes as well as crossovers between homologous chromosomes. A sexual cycle for the opportunistic human pathogenic fungus Aspergillus fumigatus is known, but the genetic consequences of meiosis have remained unknown. Among other Aspergilli, it is known that A. flavus has a moderately high recombination rate with an average of 4.2 crossovers per chromosome pair, whereas A. nidulans has in contrast a higher rate with 9.3 crossovers per chromosome pair. Here, we show in a cross between A. fumigatus strains that they produce an average of 29.9 crossovers per chromosome pair and large variation in total map length across additional strain crosses. This rate of crossovers per chromosome is more than twice that seen for any known organism, which we discuss in relation to other genetic model systems. We validate this high rate of crossovers through mapping of resistance to the laboratory antifungal acriflavine by using standing variation in an undescribed ABC efflux transporter. We then demonstrate that this rate of crossovers is sufficient to produce one of the common multidrug resistant haplotypes found in the cyp51A gene (TR34/L98H) in crosses among parents harboring either of 2 nearby genetic variants, possibly explaining the early spread of such haplotypes. Our results suggest that genomic studies in this species should reassess common assumptions about linkage between genetic regions. The finding of an unparalleled crossover rate in A. fumigatus provides opportunities to understand why these rates are not generally higher in other eukaryotes.

The long-term effects of genomic selection: 2. Changes in allele frequencies of causal loci and new mutations.

Genetic selection has been applied for many generations in animal, plant, and experimental populations. Selection changes the allelic architecture of traits to create genetic gain. It remains unknown whether the changes in allelic architecture are different for the recently introduced technique of genomic selection compared to traditional selection methods and whether they depend on the genetic architectures of traits. Here, we investigate the allele frequency changes of old and new causal loci under 50 generations of phenotypic, pedigree, and genomic selection, for a trait controlled by either additive, additive and dominance, or additive, dominance, and epistatic effects. Genomic selection resulted in slightly larger and faster changes in allele frequencies of causal loci than pedigree selection. For each locus, allele frequency change per generation was not only influenced by its statistical additive effect but also to a large extent by the linkage phase with other loci and its allele frequency. Selection fixed a large number of loci, and 5 times more unfavorable alleles became fixed with genomic and pedigree selection than with phenotypic selection. For pedigree selection, this was mainly a result of increased genetic drift, while genetic hitchhiking had a larger effect on genomic selection. When epistasis was present, the average allele frequency change was smaller (∼15% lower), and a lower number of loci became fixed for all selection methods. We conclude that for long-term genetic improvement using genomic selection, it is important to consider hitchhiking and to limit the loss of favorable alleles.

A genetic linkage map and improved genome assembly of the termite symbiont Termitomyces cryptogamus.

BACKGROUND: The termite-fungus symbiosis is an ancient stable mutualism of two partners that reproduce and disperse independently. With the founding of each termite colony the symbiotic association must be re-established with a new fungus partner. Complementarity in the ability to break down plant substrate may help to stabilize this symbiosis despite horizontal symbiont transmission. An alternative, non-exclusive, hypothesis is that a reduced rate of evolution may contribute to stabilize the symbiosis, the so-called Red King Effect. METHODS: To explore this concept, we produced the first linkage map of a species of Termitomyces, using genotyping by sequencing (GBS) of 88 homokaryotic offspring. We constructed a highly contiguous genome assembly using PacBio data and a de-novo evidence-based annotation. This improved genome assembly and linkage map allowed for examination of the recombination landscape and its potential effect on the mutualistic lifestyle. RESULTS: Our linkage map resulted in a genome-wide recombination rate of 22 cM/Mb, lower than that of other related fungi. However, the total map length of 1370 cM was similar to that of other related fungi. CONCLUSIONS: The apparently decreased rate of recombination is primarily due to genome expansion of islands of gene-poor repetitive sequences. This study highlights the importance of inclusion of genomic context in cross-species comparisons of recombination rate.

Species dynamics in natural bacterial communities over multiple rounds of propagation.

Our experimental work illustrates how microbial ecosystems can be shaped by selective pressures over long-term ecological time scales. Natural microbial ecosystems generally consist of various co-existing species, where community composition describes the frequency at which species or types are present. Overall functionality of the system is achieved by interacting species. Upon short-term selection, for instance by transfer to a novel environment, community composition and functionality may change in a process referred to as species sorting. Various factors, such as initial community composition and selective pressures from the environment, may influence this change. Mabisi is a traditional fermented food from Zambia that naturally contains a bacterial community of around twenty unique bacterial types. We used six comparable but different natural bacterial Mabisi communities, each split into five identical replicates, for 16 propagation cycles in a novel, common laboratory environment. Composition of the bacterial communities changed upon propagation. The influence of four main factors on community composition, i.e. initial composition (history), impact of the environment (adaptation), changes due to interaction between species and random processes (chance) in species dynamics, was tested using maximum likelihood ratios. Initial community composition seemed to determine the change in community composition, followed by random processes. Interestingly, we observed convergence at the level of ecosystem functionality, which was measured as profiles of metabolic output. This suggests different combinations of species or types can achieve similar eco-system functionality.

The long-term effects of genomic selection: 1. Response to selection, additive genetic variance, and genetic architecture.

BACKGROUND: Genomic selection has revolutionized genetic improvement in animals and plants, but little is known about its long-term effects. Here, we investigated the long-term effects of genomic selection on response to selection, genetic variance, and the genetic architecture of traits using stochastic simulations. We defined the genetic architecture as the set of causal loci underlying each trait, their allele frequencies, and their statistical additive effects. We simulated a livestock population under 50 generations of phenotypic, pedigree, or genomic selection for a single trait, controlled by either only additive, additive and dominance, or additive, dominance, and epistatic effects. The simulated epistasis was based on yeast data. RESULTS: Short-term response was always greatest with genomic selection, while response after 50 generations was greater with phenotypic selection than with genomic selection when epistasis was present, and was always greater than with pedigree selection. This was mainly because loss of genetic variance and of segregating loci was much greater with genomic and pedigree selection than with phenotypic selection. Compared to pedigree selection, selection response was always greater with genomic selection. Pedigree and genomic selection lost a similar amount of genetic variance after 50 generations of selection, but genomic selection maintained more segregating loci, which on average had lower minor allele frequencies than with pedigree selection. Based on this result, genomic selection is expected to better maintain genetic gain after 50 generations than pedigree selection. The amount of change in the genetic architecture of traits was considerable across generations and was similar for genomic and pedigree selection, but slightly less for phenotypic selection. Presence of epistasis resulted in smaller changes in allele frequencies and less fixation of causal loci, but resulted in substantial changes in statistical additive effects across generations. CONCLUSIONS: Our results show that genomic selection outperforms pedigree selection in terms of long-term genetic gain, but results in a similar reduction of genetic variance. The genetic architecture of traits changed considerably across generations, especially under selection and when non-additive effects were present. In conclusion, non-additive effects had a substantial impact on the accuracy of selection and long-term response to selection, especially when selection was accurate.

Somatic deficiency causes reproductive parasitism in a fungus.

Some multicellular organisms can fuse because mergers potentially provide mutual benefits. However, experimental evolution in the fungus Neurospora crassa has demonstrated that free fusion of mycelia favours cheater lineages, but the mechanism and evolutionary dynamics of this exploitation are unknown. Here we show, paradoxically, that all convergently evolved cheater lineages have similar fusion deficiencies. These mutants are unable to initiate fusion but retain access to wild-type mycelia that fuse with them. This asymmetry reduces cheater-mutant contributions to somatic substrate-bound hyphal networks, but increases representation of their nuclei in the aerial reproductive hyphae. Cheaters only benefit when relatively rare and likely impose genetic load reminiscent of germline senescence. We show that the consequences of somatic fusion can be unequally distributed among fusion partners, with the passive non-fusing partner profiting more. We discuss how our findings may relate to the extensive variation in fusion frequency of fungi found in nature.

Parasexual recombination enables Aspergillus fumigatus to persist in cystic fibrosis.

Aspergillus fumigatus is a saprobic fungus that causes a range of pulmonary diseases, some of which are characterised by fungal persistence such as is observed in cystic fibrosis (CF) patients. Creation of genetic variation is critical for A. fumigatus to adapt to the lung environment, but biofilm formation, especially in CF patients, may preclude mutational supply in A. fumigatus due to its confinement to the hyphal morphotype. We tested our hypothesis that genetic variation is created through parasexual recombination in chronic biofilms by phenotypic and genetic analysis of A. fumigatus isolates cultured from different origins. As diploids are the hallmark of parasex, we screened 799 A. fumigatus isolates obtained from patients with CF, chronic pulmonary lung disease and acute invasive aspergillosis, and from the environment for spore size. Benomyl sensitivity, nuclear content measurements through fluorescence-activated cell sorting and scanning electron microscopy were used to confirm the diploid state of large size spores. Whole genome sequencing was used to characterise diploid-associated genetic variation. We identified 11 diploids in isolates recovered from six of 11 (55%) CF patients and from one of 24 (4%) chronic aspergillosis patients, but not in 368 isolates from patients with acute Aspergillus infection and the environment. Diploid formation was associated with accumulation of mutations and variable haploid offspring including a voriconazole-resistant isolate. Parasexual recombination allows A. fumigatus to adapt and persist in CF patients, and plays a role in azole resistance development. Our findings are highly significant for understanding the genetics and biology of A. fumigatus in the human lung.

Bracon brevicornis Genome Showcases the Potential of Linked-Read Sequencing in Identifying a Putative Complementary Sex Determiner Gene.

Bracon brevicornis is an ectoparasitoid of a wide range of larval-stage Lepidopterans, including several pests of important crops, such as the corn borer, Ostrinia nubilalis. It is also one of the earliest documented cases of complementary sex determination in Hymenoptera. Here, we present the linked-read-based genome of B. brevicornis, complete with an ab initio-derived annotation and protein comparisons with fellow braconids, Fopius arisanus and Diachasma alloeum. We demonstrate the potential of linked-read assemblies in exploring regions of heterozygosity and search for structural and homology-derived evidence of the complementary sex determiner gene (csd).

The Medical Triazole Voriconazole Can Select for Tandem Repeat Variations in Azole-Resistant Aspergillus Fumigatus Harboring TR34/L98H Via Asexual Reproduction.

Azole-resistant Aspergillus fumigatus isolates recovered at high frequency from patients, harbor mutations that are associated with variation of promoter length in the cyp51A gene. Following the discovery of the TR34/L98H genotype, new variations in tandem repeat (TR) length and number of repeats were identified, as well as additional single nucleotide polymorphisms (SNPs) in the cyp51A gene, indicating that the diversity of resistance mutations in A. fumigatus is likely to continue to increase. Investigating the development routes of TR variants is critical to be able to design preventive interventions. In this study, we tested the potential effects of azole exposure on the selection of TR variations, while allowing haploid A. fumigatus to undergo asexual reproduction. Through experimental evolution involving voriconazole (VOR) exposure, an isolate harboring TR343/L98H evolved from a clinical TR34/L98H ancestor isolate, confirmed by whole genome sequencing. TR343/L98H was associated with increased cyp51A expression and high VOR and posaconazole MICs, although additional acquired SNPs could also have contributed to the highly azole-resistant phenotype. Exposure to medical azoles was found to select for TR343, thus supporting the possibility of in-host selection of TR34 variants.

Composition and Diversity of Natural Bacterial Communities in Mabisi, a Traditionally Fermented Milk.

Many traditionally fermented milk products such as mabisi involve spontaneous fermentation, which can result in bacterial community composition variation due to selection pressure. The aim of this study was to determine the composition of bacterial communities in the different types of mabisi produced across Zambia and identify the factors that influence their composition. Samples of mabisi were collected across the country, and analyzed for pH and bacterial communities using 16S rRNA amplicon sequencing. We found that the bacterial community composition was dominated by members of two phyla, i.e., Firmicutes and Proteobacteria, from which the top 10 most abundant genera were Lactococcus, Lactobacillus, Streptococcus, Enterobacter, Citrobacter, Klebsiella, Kluyvera, Buttiauxella, Aeromonas, and Acinetobacter. The most dominant genus was Lactococcus, which was present in all types of mabisi produced from all regions. The mabisi products from traditional mabisi production regions (TMPRs) were dominated by lactic acid bacteria (LAB) whereas products from non-TMPRs were dominated by non-LAB species. Tonga mabisi, the most popular type of mabisi produced in non-TMPRs, had the most complex and diverse bacterial community composition compared to the other types, which included barotse, backslopping, creamy, and thick-tonga mabisi. Other factors that influenced bacterial community composition were geographical location, fermentation duration and pH while the type of fermentation container and producer did not. This study provides new insights that can be applied in starter culture development as well as microbial functionality studies.

Genomics of sex allocation in the parasitoid wasp Nasonia vitripennis.

BACKGROUND: Whilst adaptive facultative sex allocation has been widely studied at the phenotypic level across a broad range of organisms, we still know remarkably little about its genetic architecture. Here, we explore the genome-wide basis of sex ratio variation in the parasitoid wasp Nasonia vitripennis, perhaps the best studied organism in terms of sex allocation, and well known for its response to local mate competition. RESULTS: We performed a genome-wide association study (GWAS) for single foundress sex ratios using iso-female lines derived from the recently developed outbred N. vitripennis laboratory strain HVRx. The iso-female lines capture a sample of the genetic variation in HVRx and we present them as the first iteration of the Nasonia vitripennis Genome Reference Panel (NVGRP 1.0). This panel provides an assessment of the standing genetic variation for sex ratio in the study population. Using the NVGRP, we discovered a cluster of 18 linked SNPs, encompassing 9 annotated loci associated with sex ratio variation. Furthermore, we found evidence that sex ratio has a shared genetic basis with clutch size on three different chromosomes. CONCLUSIONS: Our approach provides a thorough description of the quantitative genetic basis of sex ratio variation in Nasonia at the genome level and reveals a number of inter-related candidate loci underlying sex allocation regulation.

Nutritional Composition and Microbial Communities of Two Non-alcoholic Traditional Fermented Beverages from Zambia: A Study of Mabisi and Munkoyo.

Traditional fermented foods and beverages are common in many countries, including Zambia. While the general (nutritional) benefits of fermented foods are widely recognised, the nutritional composition of most traditional fermented foods is unknown. Furthermore, fermentation is known to add nutritional value to raw materials, mainly by adding B-vitamins and removing anti-nutritional factors. In the case of traditional fermentation, the composition of microbial communities responsible for fermentation varies from producer to producer and this may also be true for the nutritional composition. Here, we characterized the nutrient profile and microbial community composition of two traditional fermented foods: milk-based Mabisi and cereal-based Munkoyo. We found that the two products are different with respect to their nutritional parameters and their microbial compositions. Mabisi was found to have higher nutritional values for crude protein, fat, and carbohydrates than Munkoyo. The microbial community composition was also different for the two products, while both communities were dominated by lactic acid bacteria. Our analyses showed that variations in nutritional composition, defined as the amount of consumption that would contribute to the estimated average requirement (EAR), might be explained by variations in microbial community composition. Consumption of Mabisi appeared to contribute more than Munkoyo to the EAR and its inclusion in food-based recommendations is warranted. Our results show the potential of traditional fermented foods such as Mabisi and Munkoyo to add value to current diets and suggests that variations in microbial composition between specific product samples can result in variations in nutritional composition.

Distinct genomic signals of lifespan and life history evolution in response to postponed reproduction and larval diet in Drosophila.

Reproduction and diet are two major factors controlling the physiology of aging and life history, but how they interact to affect the evolution of longevity is unknown. Moreover, although studies of large-effect mutants suggest an important role of nutrient sensing pathways in regulating aging, the genetic basis of evolutionary changes in lifespan remains poorly understood. To address these questions, we analyzed the genomes of experimentally evolved Drosophila melanogaster populations subjected to a factorial combination of two selection regimes: reproductive age (early versus postponed), and diet during the larval stage ("low," "control," "high"), resulting in six treatment combinations with four replicate populations each. Selection on reproductive age consistently affected lifespan, with flies from the postponed reproduction regime having evolved a longer lifespan. In contrast, larval diet affected lifespan only in early-reproducing populations: flies adapted to the "low" diet lived longer than those adapted to control diet. Here, we find genomic evidence for strong independent evolutionary responses to either selection regime, as well as loci that diverged in response to both regimes, thus representing genomic interactions between the two. Overall, we find that the genomic basis of longevity is largely independent of dietary adaptation. Differentiated loci were not enriched for "canonical" longevity genes, suggesting that naturally occurring genic targets of selection for longevity differ qualitatively from variants found in mutant screens. Comparing our candidate loci to those from other "evolve and resequence" studies of longevity demonstrated significant overlap among independent experiments. This suggests that the evolution of longevity, despite its presumed complex and polygenic nature, might be to some extent convergent and predictable.

Fermented cereal-based Munkoyo beverage: Processing practices, microbial diversity and aroma compounds.

Fermented cereal-based foods play a crucial role in attaining food and nutrition security for resource-poor populations in sub-Saharan Africa. These products are widely produced by spontaneous fermentation using of cereal grains as raw material. They have a unique taste and flavour, are rich sources of energy and their non-alcoholic nature makes them ideal for consumption by the entire population, including children. Lactic acid bacteria dominate the fermentation process and lead to a low pH of around 4, which suppresses the growth of pathogenic bacteria, thereby increasing the shelf-life and safety of the food. Knowledge about processing practices, consumption patterns and bacterial communities is essential to regulate processing and design appropriate mixes of micro-organisms to produce starter cultures for commercial production of standard-quality fermented foods that meet desired quality characteristics. In four regions of Zambia, we surveyed processing practices and consumption patterns of a spontaneously fermented cereal-based beverage called Munkoyo, commonly produced in Zambia and the Democratic Republic of Congo. Variations in processing practices exist in cooking time of the unfermented maize porridge and time allowed for fermentation. Consumption is mainly at household level and the product is considered as an energy drink. Characterisation of the bacterial communities of over 90 samples with 16S amplicon sequencing on DNA extracted from the entire bacterial community revealed six dominant families, namely Streptococcaceae, Leuconostocaceae, Enterobacteriaceae, Lactabacillales, Bacillaceae and Aeromonadaceae, and a Shannon index of up to 1.18 with an effective number of 3.44 bacterial species. Bacterial communities that underlie the fermentation in Munkoyo differ in their composition for the different regions using common processing steps, suggesting that different combinations of bacteria can be used to achieve successful Munkoyo fermentation. Analysis of aroma profiles in 15 different samples from two different Provinces showed that aldehydes, esters, organic acids, alkanes, alkenes and alcohols dominated.

The genome of the live-bearing fish Heterandria formosa implicates a role of conserved vertebrate genes in the evolution of placental fish.

BACKGROUND: The evolution of complex organs is thought to occur via a stepwise process, each subsequent step increasing the organ's complexity by a tiny amount. This evolutionary process can be studied by comparing closely related species that vary in the presence or absence of their organs. This is the case for the placenta in the live-bearing fish family Poeciliidae, as members of this family vary markedly in their ability to supply nutrients to their offspring via a placenta. Here, we investigate the genomic basis underlying this phenotypic variation in Heterandria formosa, a poeciliid fish with a highly complex placenta. We compare this genome to three published reference genomes of non-placental poeciliid fish to gain insight in which genes may have played a role in the evolution of the placenta in the Poeciliidae. RESULTS: We sequenced the genome of H. formosa, providing the first whole genome sequence for a placental poeciliid. We looked for signatures of adaptive evolution by comparing its gene sequences to those of three non-placental live-bearing relatives. Using comparative evolutionary analyses, we found 17 genes that were positively selected exclusively in H. formosa, as well as five gene duplications exclusive to H. formosa. Eight of the genes evolving under positive selection in H. formosa have a placental function in mammals, most notably endometrial tissue remodelling or endometrial cell proliferation. CONCLUSIONS: Our results show that a substantial portion of positively selected genes have a function that correlates well with the morphological changes that form the placenta of H. formosa, compared to the corresponding tissue in non-placental poeciliids. These functions are mainly endometrial tissue remodelling and endometrial cell proliferation. Therefore, we hypothesize that natural selection acting on genes involved in these functions plays a key role in the evolution of the placenta in H. formosa.

Adaptation to developmental diet influences the response to selection on age at reproduction in the fruit fly.

Experimental evolution (EE) is a powerful tool for addressing how environmental factors influence life-history evolution. While in nature different selection pressures experienced across the lifespan shape life histories, EE studies typically apply selection pressures one at a time. Here, we assess the consequences of adaptation to three different developmental diets in combination with classical selection for early or late reproduction in the fruit fly Drosophila melanogaster. We find that the response to each selection pressure is similar to that observed when they are applied independently, but the overall magnitude of the response depends on the selection regime experienced in the other life stage. For example, adaptation to increased age at reproduction increased lifespan across all diets; however, the extent of the increase was dependent on the dietary selection regime. Similarly, adaptation to a lower calorie developmental diet led to faster development and decreased adult weight, but the magnitude of the response was dependent on the age-at-reproduction selection regime. Given that multiple selection pressures are prevalent in nature, our findings suggest that trade-offs should be considered not only among traits within an organism, but also among adaptive responses to different-sometimes conflicting-selection pressures, including across life stages.

Selective Survival of Embryos Can Explain DNA Methylation Signatures of Adverse Prenatal Environments.

An adverse intrauterine environment is associated with long-term physiological changes in offspring. These are believed to be mediated by epigenomic marks, including DNA methylation (DNAm). Changes in DNAm are often interpreted as damage or plastic responses of the embryo. Here, we propose that stochastic DNAm variation, generated during remodeling of the epigenome after fertilization, contributes to DNAm signatures of prenatal adversity through differential survival of embryos. Using a mathematical model of re-methylation in the early embryo, we demonstrate that selection, but not plasticity, will generate a characteristic reduction in DNAm variance at loci that contribute to survival. Such a reduction in DNAm variance was apparent in a human cohort prenatally exposed to the Dutch famine, illustrating that it is possible to detect a signature of selection on epigenomic variation. Selection should be considered as a possible mechanism linking prenatal adversity to subsequent health and may have implications when evaluating interventions.

The obligate alkalophilic soda-lake fungus Sodiomyces alkalinus has shifted to a protein diet.

Sodiomyces alkalinus is one of the very few alkalophilic fungi, adapted to grow optimally at high pH. It is widely distributed at the plant-deprived edges of extremely alkaline lakes and locally abundant. We sequenced the genome of S. alkalinus and reconstructed evolution of catabolic enzymes, using a phylogenomic comparison. We found that the genome of S. alkalinus is larger, but its predicted proteome is smaller and heavily depleted of both plant-degrading enzymes and proteinases, when compared to its closest plant-pathogenic relatives. Interestingly, despite overall losses, S. alkalinus has retained many proteinases families and acquired bacterial cell wall-degrading enzymes, some of them via horizontal gene transfer from bacteria. This fungus has very potent proteolytic activity at high pH values, but slowly induced low activity of cellulases and hemicellulases. Our experimental and in silico data suggest that plant biomass, a common food source for most fungi, is not a preferred substrate for S. alkalinus in its natural environment. We conclude that the fungus has abandoned the ancestral plant-based diet and has become specialized in a more protein-rich food, abundantly available in soda lakes in the form of prokaryotes and small crustaceans.